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The collective behavior of a coupled map lattice having unbounded chaotic local dynamics is investigated through the properties 
of its mean field. The presence of unstable periodic orbits in the local maps determines the emergence of nontrivial collective 
behavior. Windows of collective period-two states are found in parameter space. 



PACS numbers: 05.45. +b, 02.50.- 



Coupled map lattices (CML) have provided fruitful 
models for the study of many spatiotemporal processes in 
a variety of contexts Recently, there has been great 
interest in the use of CML models in the investigation of 
global behaviors of networks of chaotic elements [||,|],f§ ■ 
Such collective phenomena have many implications rang- 
ing from the fundamentals of statistical mechanics |^| to 
biological information processing, and even possible prac- 
tical applications ||. 

In this letter, we study the collective behavior of the 
one-dimensional CML system 



x t+l {i) = f{x t {i)) + j[f(x t (i - 1)) + f(x t (i + 1)) 
-2/(x t (*))], 



(1) 



where t is a discrete time, i labels the lattice sites 
(i = 1 , . . . , , N) , 7 is a parameter measuring the diffusive 
coupling strength between neighboring sites, and / is a 
map describing the local dynamics. Periodic boundary 
conditions are assumed. 

The collective behavior of the lattice can be monitored 
through the instantaneous mean field or activity of the 
network, defined as the space average of the local vari- 
ables x(i) at time t M: 



1 N 



(2) 



A wide variety of local map functions / have been em- 
ployed in CML models depending on the particular appli- 
cation. Usually, bounded maps belonging to some uni- 
versality class are considered to be the source of local 
chaos in the study of collective spatiotemporal dynamics 
in CMLs with local |^|| or global couplings [Q. Here 
we investigate the collective behavior of the lattice de- 
scribed by (|lj) having unbounded local dynamics such as 
the logarithmic map 



f(x) = 6 + ln|x|. 



(3) 



This map posseses no maximum or minimum and its 
Schwarzian derivative is always positive. Two stable 
fixed points satisfying f(x*) = x* exist for this map: 



x\ < — 1, for b < — 1, which becomes unstable at b = — 1; 
and xZ, > 1, for b > 1, which arises from a tangent bifur- 
cation at b = 1. Chaos occurs in the parameter region 
b G [—1,1]. There are no gaps separating chaotic bands 
at any given value of & 6 [—1,1] and no periodic windows 
appear in any subinterval of b in this range §. These 
special features, contrasting with the universal proper- 
ties of the commonly used local maps, make maps of the 
type (||) an alternative tool for exploring the essential 
requirements for the emergence of collective behavior in 
CMLs and the generality of this phenomenon. 

The dynamics of the CML (Q) can be written in vector 
form as 



x t+1 = f(x t )+ 7 Mf(x 4 ), 



(4) 



where the N-dimensional vectors x t and f(x t ) have com- 
ponents [x t ]i = x t (i) and [f(x t )]i = f(x t (i)), respec- 
tively; and M is a N x N tridiagonal, symmetric ma- 
trix expressing diffusive coupling among the components 
[x.t]i- The non-zero components of M are Ma = —2; 
My = 1 {i=j±l). 

With the local map (0), the stable, spatially homo- 
geneous, stationary states x t (i) — x*, for b < — 1; or 
x t (i) = x%, for b > 1, Vi, are possible for the system (Q). 
The linear stability analysis of these states leads to the 
bifurcation conditions 



(l + 7Mfc)/'« 2 ) = ±l, 



(5) 



where {^k ■ k = 1, . . . , N} is the set of eigenvalues of 
the coupling matrix M which satisfy fik S [—4,0] ||. 
Equations (^) yield boundary curves in the parameter 
plane (7, b) which determine where each homogeneous, 
stationary state can be observed. Figure 1 shows the 
first stability boundaries for each state, corresponding to 
fik — —4. The state Xt(i) — x*, Vi, is stable for pa- 
rameter values in the region enclosed by the r.h.s. ±1 
boundaries corresponding to x\ in Eq. (§), with b < -1. 



Similarly, the state x t (i) 
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Vi, is stable inside the 



parameter region bounded by the r.h.s. ±1 curves asso- 
ciated to x\, with b > 1. Within these regions of stable 
homogeneous stationary states, the asymptotic values of 
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the mean field are St = x* and St — x\ , respectively. The 
crossing of either boundary signals the appearance of a 
spatially inhomogeneous state which should be reflected 
in a dispersion of St- 

When the value of the parameter b is in the range 
[—1,1], corresponding to local chaotic dynamics, the 
asymptotic collective behavior of the lattice, as given 
by St->oc , reveals the existence of global low-dimensional 
periodic attractors, subjected to fluctuations of intrinsic 
statistical origin. Figure 2 shows the bifurcation dia- 
gram of the asymptotic mean field St as a function of 
the local parameter b. The quantity St was calculated 
at each time step during a run starting from random ini- 
tial conditions uniformly distributed on the interval [-8,4] 
at each site and for each value of b, after discarding the 
transients. The coupling strength was fixed at 7 = 0.5, 
corresponding to a "totalistic" coupling of the "game of 
life" type 0. Similar bifurcation diagrams can be ob- 
tained for other fixed values of the coupling parameter. 

For b < — 1 or & > 1, the asymptotic mean field is iden- 
tical to the values of the corresponding fixed points of the 
single map (|^) in these ranges of b, as expected from the 
stability diagram of the homogeneous stationary states 
in Fig. 1. In the region b G [—1,1], figure 2 shows a 
pitchfork bifurcation at a value b c ~ —0.52 from a col- 
lective fixed point state (a state for which the time series 
of St manifestly displays statistical fluctuations around 
a single value) to a collective period-two state (a state 
for which the time series of St alternatively varies be- 
tween the corresponding neighborhoods of two separated 
values). The small vertical segments seen in Fig. 2 can 
be interpreted as the amplitude of the fluctuations of 
the mean field about the corresponding global attractor 
at given parameter values. The fluctuations around the 
global stationary or the period-two attractors are due to 
the fact that the local variables behave chaotically, as 
can be attested by the time series of any site and by the 
existence of a positive largest Lyapunov exponent. 

The appearance of a period-two collective state in the 
one-dimensional lattice (0) is related to the fact that the 
iterates of the local map (|3|) move alternatively from val- 
ues above the unstable fixed point x\ to values bellow 
this point in the interval b £ [— 1, b c ], even though there 
are no separated chaotic bands 0. The iterates behave 
chaotically within each side of the unstable fixed point 
x\, but periodically about it. The unstable fixed point x\ 
establishes a "symmetry" line around which the iterates 
oscillate with period two. As a comparison, such effect 
does not occur in the unbounded map g(x) = a — 1/x, 
which lacks unstable periodic orbits for a £ [—2, 2]; con- 
sequently, no collective periodic behavior emerges in a 
coupled lattice of these maps at that range of the param- 
eter a ||. 

Coupling induces synchronization of the array in the 
range b e [— 1,6 C ], in the sense that iterates of the lo- 
cal chaotic sites tend to move together above and bellow 



the unstable fixed point x\ at alternate time steps. The 
individual local values may be different within each of 
these two regions at a given instant. Figure 3(a) presents 
the asymptotic mean field St as a function of the cou- 
pling strength 7, keeping constant b = —0.8. When 7 is 
close to zero, and starting each time from random ini- 
tial conditions, St fluctuates about an average value of 
the desynchronized local maps. There exists a critical 
value of the coupling at 7 C ~ 0.01 where the transition 
from a statistical average to a period-two collective state 
takes place. The periodic behavior of St remains practi- 
cally unchanged until the coupling reaches a second crit- 
ical value Yc — 0-53 at which synchronization is again 
lost. Figure 3(b) shows the asymptotic behavior of an 
arbitrary individual site simultaneously monitored as a 
function of the coupling. The value of the unstable fixed 
point x\ corresponding to b = —0.8 is indicated by a 
horizontal line. Coupling enhances the separation of al- 
ternate iterates about x\, creating a gap between them 
for 7 G [7c, 7c]- The mean field reflects this separation 
as well as the synchronization of the lattice produced by 
the coupling. 

Figure 4 shows the dependence of the asymptotic St 
on the size of the lattice N, with the other parameters 
held fixed (7 = 0.5, b = —0.8). There is a rather small 
critical size N c ~ 10 at which the global period-two at- 
tractor distinctively emerges. As in Fig 2, the vertical 
segments represent fluctuations around the global period- 
two attractor and they correspond to even or odd steps 
of the asymptotic time series of St, respectively. The 
variance of either subset of steps of the time series of St 
decreases as A^ _1 , as expected. However, the variance 
of the mean field itself, for the chosen parameters, tends 
to a constant value for large enough N. In the limit 
N — > 00, the global period- two orbit will have the form 

St ■ ■ ■ ■ Si, S2, Si, S2 The time average of the mean 

field will be S — (si + S2)/2, while the variance a will 
yield 

t=i t=i 

For the parameters of Fig. 4, s\ ~ —2.55, S2 — 0.08. The 
limiting value a ~ 1.72 is approached for N ~ 10 3 . The 
existence of a saturation value for the variance of the 
mean field characterizes the emergence of nontrivial (i.e. 
periodic, quasiperiodic, or chaotic) collective behavior. 
This phenomenon has been called "violation of the law 
of large numbers" in the context of globally coupled maps 
belonging to some universality class (tent, quadratic, or 
circle maps) [flffo] . 

The onset of a periodic collective state at some values 
of the parameters of the system is reminiscent of a phase 
transition, because it corresponds to abrupt changes in 
the statistical properties of the lattice, as described by 
St->oc, a quantity which acts as an order parameter. We 
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have found windows of global periodic behavior in the 
space of parameters. Periodic collective states have been 
observed in higher dimensional lattices of locally chaotic 
coupled maps. In those cases, the collective behavior con- 
sists in statistical cycling of the mean field (or the prob- 
ability density) among chaotic bands of the local one- 
hump maps fllQ. This paper shows that high space 
dimension, large system size, strong coupling, bounded 
iterates, gap-separated chaotic bands, or the existence of 
periodic windows in the local dynamics, are not essential 
requirements for the emergence of nontrivial collective 
behavior. We have found that coupling can synchronize 
the cycling of the chaotic iterates and enhances their sep- 
aration around unstable periodic orbits of the local maps 
in some parameter ranges giving rise to periodic global 
behavior in a lattice. Our results suggest that the emer- 
gence of nontrivial collective behavior should be a rather 
generic phenomenon in deterministic systems of coupled 
chaotic elements, where unstable periodic orbits are al- 
ways present. 
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FIG. 1. Stability boundaries for the homogeneous station- 
ary states of the CML (Q). The labels on each curve identify 
the first stability boundaries with ±1 in Eq. (||) for each state, 
corresponding to u,k — —4. The vertical lines at b = — 1 and 
6=1 correspond to u,u = 0. Curves associated with other 
eigenvalues lie outside the regions enclosed by the boundaries 
shown. The states St = x\ and St = x\ are stable inside the 
indicated regions. 

FIG. 2. Bifurcation diagram of the asymptotic St, as a 
function of the local parameter b. Coupling is fixed at 7 = 0.5. 
For each value of b, 100 iterates are shown, after discarding 
5000 transients. Lattice size is N = 10 5 . For b < -1 and 
b > 1, St corresponds to the values of the homogeneous sta- 
tionary states x* and a;|, respectively. 

FIG. 3. a) St as a function of the coupling parameter 7. 
Local parameter is fixed at b — —0.8. For each value of 7, 100 
iterates are shown, after discarding 5000 transients. Lattice 
size is TV = 10 5 . b) Asymptotic iterates of site i — 100, for the 
same parameters as in a). The horizontal line corresponds to 
the value of the unstable x\ for the given b. 

FIG. 4. Asymptotic mean field St as a function of lattice 
size N. Fixed parameters are 7 — 0.5; b — —0.8. 
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